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PRELIMINARY  PARAMETRIC  STUDY 
OF  THE  FLUTTER  ARMING  PRINCIPLE 
FOR  FUZES 

Prepared  by: 

P.  D.  Grattpn 


ABSTRACT:  The  Flutter  Arming  Mechanism  (FAM)  principle  is  a 
controlled  flutter  phenomenon.  A  rectangular  flat  plate  oscillating 
member  is  pivoted  about  its  midchord  on  a  member  providing  a  flexural 
restoring  moment.  Placed  edgewise  in  an  air  stream,  the  system  is  in 
unstable-stable  equilibrium;  at  a  predetermined  air  speed  or  above, 
aerodynamic  lift  on  the  front  portion  of  the  flat  plate  overcomes  the 
flexural  resisting  moment  and  the  oscillator  vibrates.  The  vibrations 
occur  at  the  natural  frequency  of  the  oscillator-flexural  member 
combination.  Wind  tunnel  data  are  presented  Which  indicate  the  FAM 
principle  has  the  necessary  basic  criterion  of  a  time  base,  velocity 
discriminating  mechanism  required  for  bomb  and  bomblet  fuze 
application. 


Published 

11  SEPTEMBER  1967 


U.  S.  NAVAL  ORDNANCE  LABORATORY 
WHITE  OAK,  MARYLAND 


i 


NOLTR  67-144  11  September  1967 

PRELIMINARY  PARAMETRIC  STUDY  OF  THE  FLUTTER  ARMING  PRINCIPLE  FOR  FUZES 

This  report  is  the  first  in  a  series  of  reports  to  be  generated 
covering  the  parameters  and  functioning  of  the  flutter  arming 
principle.  It  represents  the  first  parametric  study  conducted  with 
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INTRODUCTION 

1.  The  components  m  a  conventional  mechanical  fuzing  system 
for  bombs  include  an  arming  vane*  to  extract  energy  from  the  envi¬ 
ronment,  high  speed  bearings  to  insure  proper  rotational  character¬ 
istics  of  the  arming  vanes,  clutch  or  similar  mechanism  to  sense 

a  specific  velocity  threshold  and  a  clock  escapement  or  governor  for 
proper  delay  arming  of  the  fuze.  This  technical  report  is  concerned 
with  another  mechanism  that  current  exploratory  development  shows  to 
have  promise  for  accomplishing  the  same  ends.  At  a  minimum  velocity 
and  above,  controlled  vibrations  from  a  flat  plate  pivoted  about  its 
midchord  provide  3  constant  frequency  time  base  mechanism.  The 
arming  vanes,  hi^h  speed  bearings  and  clock  escapement  or  governor 
are  eliminated.  The  mechamism  has  the  potential  of  becoming  an 
inexpensive  environmentally  operated  delay  arming  device  that  could 
effect  a  considerable  saving  when  used  on  fin  stabilized  ordnance. 

2.  This  fuze  arming  mechanism,  which  is  based  on  the  phenomenon 
of  aerodynamic  flutter,  is  known  as  the  Flutter  Arming  Mechanism 
(FAM) .  Studies  of  the  FAM  were  conducted  in  the  wind  tunnel  by  the 
use  of  a  special  model  of  the  FAM  referred  to  as  the  control  volume 
model.  This  model  was  so  constructed  that  certain  fundamental  para¬ 
meters  could  be  changed  from  test  to  test  to  study  the  effects  of 
the  change.  This  model  was  not  intended  to  represent  any  scale  fuze 
design  but  idealizes  the  operational  environment  and  air  flow 
considerations.  The  information  obtained  from  this  model  will  be 
applied  to  fuze  design  models  in  the  future.  The  majority  of  the 
data  presented  were  for  subsonic  tunnel  speeds  but  some  supersonic 
runs  were  conducted. 


SYMBOLS 

3.  For  convenience  in  reference  all  symbols  used  in  subsequent 
paragraphs  are  defined  here. 

a  =  distance  from  pivot  point  to  center  of  pressure  (in.) 
a*  =  speed  of  sound  at  throat  (ft/’sec) 
f  =  natural  frequency  (Hz) 
qc  =  chamber  dynamic  pressure  (psi) 

Ac  =  average  chamber  area  (in.) 

Ac  =  chamber  area  (in.) 

=  inlet  area  (in.) 

At  =  throat  area  (in.) 

CN  =  aerodynamic  lift  coefficient 

Ew  =  Young's  Modulus,  30  x  10^  psi  2 

I  =  oscillator  moment  of  inertia  (lb-in. -sec  ) 

Iw  =  inertia  of  wire  (in. 4) 

K  =  spring  constant  (lb-in. /rad  ) 

*  typically  a  rotating  propeller  or  turbine,  axial  or  radial  flow 
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M  =  Mach  number 

=  Aerodynamic  Moment 

Pc  =  chamber  static  pressure  at  leading  edge  of  oscillator 

S  =  surface  area  of  oscillator  (in. 2) 

Vc  =  chamber  velocity  at  leading  edge  of  oscillator  (fV'sec) 
Vi  =  inlet  or  upstream  air  velocity  (ft/sec) 

Vcr  =  chamber  velocity  at  threshold  (ft/sec) 

Z  =  length  of  wire  between  supports  (in.) 

«  =  angle  of  attack 

p_  =  chamber„air  density  at  leading  edge  of  oscillator 
C  ( lb-sec2) /ft* 

=  inlet  air  density  (lb-sec2)/ft4 
=  inlet  static  pressure  (cm  of  mercury) 

Pfci  =  inlet  total  pressure  (cm  of  mercury) 

Ptc  =  chamber  total  pressure  (cm  of  mercury) 
t  -  temperature  °  C 

d  =  location  of  oscillator  leading  edge  from  nozzle 
exit  (in.) 


THEORY 

4.  The  unstable  divergent  motion  witnessed  by  wing  and  tail 
sctions  on  some  aircraft  is  a  type  of  vibration  that  once  started 
t  a  critical  speed  the  end  result  is  an  ultimate  structural 
lilure.  Another  extreme  would  be  the  dartped  motion  where  at  first 

member  would  vibrate  when  excited  but  internal  and  external 
inping  causes  gradual  subsiding  of  motion  in  a  finite  amount  of 
Lme.  The  FAM  principle  is  an  unstable  oscillating  type  or 
>ntrolled  flutter  mecbamism  between  these  two  extremes  (continuous 
jtion  with  no  change  in  amplitude) . 

5.  The  control  volume  model  has  a  flat  plate  or  oscillator 
.voted  and  balanced  about  its  midchord  located  in  a  parallel  walled 
lamber.  A  restoring  spring  is  attached  at  the  pivot  axis  to  give 
bias  zero  point  of  reference.  An  air  stream  inpinging  the  leading 
Ige  of  the  oscillator  will  create  an  aerodynamic  lifting  force  that 
>rks  at  the  center  of  pressure  (c.p.)  tending  to  rotate  the  oscil- 
tor  from  the  zero  reference  point.  Only  at  a  threshold  speed  and 
>ove  will  the  aerodynamic  lifting  forces  be  sufficient  to  overcome 
«e  spring  restoring  moment  causing  the  oscillator  to  rotate  to  a 
all  position  (approximately  10  degrees) .  At  stall  the  aerodynamic 
ft  is  essentially  zero  but  the  restoring  spring  has  acquired  inter- 
1  energy  (due  to  bending)  that  causes  the  oscillator  to  return  to 

e  reference  point.  However,  the  oscillator  inertia  causes  a  slight 
ershoot  past  the  reference  point  allowing  the  aerodynamic  forces 
build  up  on  the  reverse  face.  The  oscillator  rotates  to  a  -10 
gree  position  (stall)  and  a  controlled  cyclic  pattern  of  amplitude 
d  frequency  is  established.  The  frequency  of  vibration  should  be 
e  natural  undamped  frequency  of  the  oscillator-restoring  member 
mbination.  A  balance  between  the  oscillator  inertia,  elastic 
operties  of  the  restoring  member,  and  the  aerodynamic  forces  can 
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give  various  threshold  velocities  for  starting  vibration  along  with 
constant  frequencies  for  a  time  base.  The  described  system  is  a 
simplified  explanation  of  the  FAM  operation.  In  actual  application 
the  control  volume  model  has  a  convergent  nozzle  that  directs  the 
air  flow  to  the  leading  edge  of  the  oscillator.  The  aerodynamic 
forces  applied  to  the  osci llator  are  chanqed  by  the  addition  of  the 
nozzle.  The  lift  is  destroyed  at  ''stall”  due  in  part  to  the 
oscillator  hiding  behind  part  of  the  nozzle  or  out  of  the  effective 
air  flow. 


SETUP 

6.  Figure  1  shows  the  relative  size  and  location  of  the  various 
components  of  the  control  volume  model.  The  housing  was  mounted  on 
a  15"  sting  and  this  in  turn  was  mounted  on  a  rotating  sector  arm 

in  the  wind  tunnel.  The  purpose  of  this  setup  was  to  rotate  the 
control  volume  model  through  angles  of  attack  yet  keep  the  leading 
edge  of  the  oscillator  in  the  center  of  the  wind  tunnel  section. 

7.  The  plane  of  the  oscillator  in  the  control  volume  model  is 
normally  horizontal;  this  is  referenced  to  a  zero  angle  of  attack. 

A  20  degree  angle  of  attack  is  characterized  when  the  control  volume 
model  is  rotated  20  degrees  from  the  horizontal.  Figure  2  shows  the 
size  and  configuration  of  the  blunt  oscillator  used  for  the  majority 
of  the  tests. 

8.  The  restoring  member  (':016  diameter  music  wire)  vas  used  in 
all  the  tests.  Wire  breakage  was  common  after  an  average  of  four  to 
five  tunnel  runs  (about  7.5  seconds  per  run),  especially  at  the 
high  subsonic  Mach  numbers. 

9.  The  channeling  or  nozzling  ’ s  needed  to  obtain  a  constant 
frequency  of  vibration  over  the  velocity  and  angle  of  attack  range. 
The  need  has  been  demonstrated  in  previous  tests  (ref.  (a)).  Various 
nozzle  openings  were  inserted  and  adjusted  to  change  the  inlet-to- 
throat  and  c’namber-to-throat  ratios. 

PARAMETERS 

10.  The  major  parameters  that  influence  the  functioning  of  the 
oscillator  are  the  area  opening  ratios  and  the  distance  from  the 
nozzle  exit  aft  to  the  leading  edge  of  the  osci llator.  The  area 
opening  ratios  refer  to  the  inlet-to-throat  areas  (A^/A^)  and  the 
average  chamber- to-throat  areas  (A^/A^) .  These  variables  were 
readily  changeable  with  slide-in-type  nozzles  (oscillator  remaining 
in  the  same  position).  For  space  locations  see  Figure  3. 

11.  The  inlet  air  speed  was  varied  from  75  to  690  knots  to 
determine  frequency  response  and  threshold  velocities. 

12.  The  angle  of  attack  on  the  oscillator  in  the  control  volume 
model  was  swept  through  angles  of  attack  starting  from  -6  and  endi ng 
at  20  degrees.  This  range  was  the  physical  limit  of  the  sector  arm 
in  the  tunnel  and  represents  an  adequate  fuze  functioning  range. 

13.  The  majority  of  tests  had  a  blunt  leading  edged  oscillator 
but  a  few  runs  were  made  with  a  rounded  leadinq  edge,  the  radius 
being  one-half  the  oscillator  thickness.  (See  Figure  2.) 
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INSTRUMENTATION 

14-  For  frequency  measurements  «  magnetic  transducer  was  inserted 
n  the  housing  just  over  the  aft  end  of  the  oscillator.  As  the  oscil- 
ator  moved  through  the  flux  lines  set  up  by  the  magnet,  a  flux 
hange  and  in  turn  voltage  change  were  permanently  recorded  by  two 
ystems.  For  Immediate  read  out,  an  X-Y  plotter  was  used.  For  more 
ophisticated  and  accurate  information  the  data  was  put  on  tape 
Digital  Analog  Recording  Equipment)  to  be  analyzed  at  a  later  date, 
he  frequency  was  recorded  as  a  function  of  time  and  angle  of  attach. 

15.  Oscillator  rise  time  was  measured  using  the  transducer  and 
visocorder.  A  carrier  frequency  0.00  Hz)  was  recorded  simultane- 

usly  with  the  transducer  output  on  light  sensitive  paper, 
omparing  the  two  outputs  gave  the  corresponding  rise  time. 

16.  Amplitude  of  oscillation  was  recorded  on  film  using  a  dial 
ndicator  attached  to  the  oscillator  shaft.  The  film  could  be  later 
nalyzed  for  exact  displacements.  An  "E"  core-type  system  will  be 
sed  in  the  future  to  record  both  amplitude  and  frequency  simultane- 
usly,  Visual  amplitude  was  constant  at  +10°. 

17.  Wind  tunnel  supply  pressures  P^,  tunnel  static  pressures 
L  and  temperatures  t  were  also  measured  and  recorded  with  Digital 
lalog  Equipment.  The  control  volume  model  chamber  static  pressure 

.  was  measured  with  a  pressure  tap  located  half  in  the  air  directly 
i  front  of  the  oscillator  and  the  other  half  of  the  tap  over  the 
scillator . 


CALCULATIONS 


18.  The  calculated  frequency  of  oscillator  vibration  is  the 
itural  frequency  of  the  oscillator  inertia  and  the  restoring  member 
nribination.  This  frequency  is  for  a  no-load  output;  no  gears, 
itchet  or  power  take-off  from  the  oscillator  to  perform  the  fuze 
tfing  and  arming  function  will  be  considered  at  this  time.  The 
equency  would  adjust  to  a  lower  value  if  the  gear  inertia,  ratchet 
ficiency  and  bearing  friction  are  considered  in  the  oscillator 
«rtia  term. 


E 

I1 

where  K 

=  12  — 

w 

Z 

for  , 

.016  dia. 

Iw  =  .322  x  10_8 

Ew  =  30  x  106  psi 

Z  =  1.140  in 

K  = 

102.0  x 

10“2lb-in 

Equation  (1) 


4 


NOLTR  67-144 


5 


NOLTR  67-144 


'(•«/*;)%  -Mi) '[‘Mil 


'  r*  -  *(s)  ] 


Equation  (5) 


which  is  in  quadratic  form  and  easily  solvable.  Attention  is  fixed 
to  the  positive  answer  only 

v  =  ~k  +Vb2  -  4ac 
c  2a 

20.  The  Ac  term  used  in  Equation  5  and  noted  on  the  Figures 
diffea  slightly  from  the  term  A  .  The  difference  is  slight  but 
for  accurate  Vc  calculations  itcmust  exist.  The  chamber  static 
pressure  probe  located  over  the  oscillator  and  partly  in  the  free 
stream  reads  an  average  pressure  over  these  two  extremes.  The 

Sc  term  is  therefore  an  average  chamber  area  corresponding  to  the 
average  static  pressure  reading. 

21.  The  chamber  stagnation  pressure  Ptc  can  be  calculated  using 
the  isentropic  equation  behind  the  shock  wave, 


p  =  p 

tc  c 


1  -  1/6  ^ 


Equation  (6) 


22.  The  inlet  and  chamber  dynamic  air  densities  were  calculated 
from  a  thermally  perfect  dry  air  equation,  (ref.(b)) 


p  =  0.-P,01293  H -  grama - Equation  (7) 

t  1  +  0.00367  t  76  milliliter^  . 


where  H  =  inlet  (P  ,)  or  chamber  (P.  )  total 

pressure^Tcm)  c 

t  =  temperature  °C 

23.  Calculating  the  Pc/Ptc  ratio  and  refering  to  ref.(c)a  Mach 
number  in  the  chamber  and  a  ratio  of  static  density/dynamic  density 
(Pc/Ptc)  can  be  obtained.  The  static  chamber  density  can  then  be 
found  by  multiplying  the  Pc/Ptc  ratio  by  Ptc  (Equation  7).  The 
static  chamber  density  was  used  in  calculating  the  threshold  or 
starting  velocity. 
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24.  The  threshold  or  start-up  velocity  can  be  calculated  knowing 
the  restoring  spring  physical  properties  and  the  aerodynamic  moment 
causing  deflection. 

Critical  or  Start-up  velocity: 

Restoring  Spring  Moment  =  Aerodynamic  Moment 

Ka  =fCN  2  PcVcr)S*  Equation  (8) 


Cjj  =  2TTa  for  a  flat  plate 


K  «  102.0  x  10“2  lb-in. 
S  =  1.19  in2 


a  =  .465  in 

Vcr  =  6,25  y  2.16  x  TCrT  .“465 

VcJ.  *  205.  ft/sec  or  121  Knots  chamber  velocity 
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25.  Certain  parameters  regarding  the  aerodynamic  lift  should  be 
mentioned. 

a.  The  lift  coefficient,  Cjj,  for  a  flat  plate  is  in  general 
2na  for  ideal  flow  conditions.  The  flow  is  not  ideal  in  this  appli¬ 
cation  due  to  the  diffuser  effect  as  the  air  enters  the  oscillator 
chamber.  The  disturbed  air  flow  along  with  the  oscillator  thickness 
and  leading  edge  configuration  will  have  some  influence  on  the  lift 
coefficient. 


b.  Another  general  assumption  that  has  been  made  deals  with 
the  center  of  pressure,  a,  location.  A  quarter  chord  from  the  lead¬ 
ing  edge  is  the  assumed  location,  regardless  of  the  angle  of  attack. 
Like  all  wing  sections  the  c.p.  for  the  flat  plate  will  shift  when 
the  angle  of  attack  changes. 

c.  The  magnitude  of  the  aerodynamic  lift  and  moment  on  the 
oscillator  can  be  altered  by  the  presence  and  location  of  the  confin¬ 
ing  chamber  walls.  Present  information  (ref .(  e))  indicatesa  maximum 
of  32%  increase  in  lift  can  be  realized  for  an  oscillator  chord  to 
chamber  wall  separation  ratio  of  1.  The  trend  is  toward  an  increas¬ 
ing  lift  coefficient  for  greater  chord/wall  ratios.  The  information 
did  not  include  the  physical  limits  of  the  control  volume  model 
configuration  of  chord/wall3 2 .8  so  this  correction  will  not  be  made. 
The  calculated  threshold  chamber  velocity  is  121  knots. 

d.  The  compressibility  factor  of  air  was  not  taken  into 
account  for  Equation  8.  The  error  amounts  to  approximately  2%  of 
the  threshold  at  121  knots  chamber  velocity  but  will  be  significantly 
greater  percentage  of  error  at  the  high  subsonic  speeds. 

26.  At  a  fixed  tunnel  air  speed  or  inlet  air  speed,  the  air  is 
accelerated  as  it  approaches  the  throat  of  the  convergent  nozzle. 

The  air  flow  is  accelerated  to  a  limiting  value  of  the  speed  of 
sound  or  Mach  1  at  the  throat.  Using  inlet  and  throat  dimensions 
for  the  nozzle,  the  inlet  air  velocity  (Vf)  establishing  Mach  1  at 
the  throat  can  be  calculated  as  shown  by  Equation  (9) . 


For  Mach  1  at  the  throat: 


v 


a* 

A±/a 


t 


Equation  (9) 


V 


i 


a*  =  1010  ft/sec 
=  =  180.5  ft/e-c 


Ai/At 


5.6 


Vi  =  107  knots 


8 


I 


NOLTR  67-144 
RESULTS 

27.  Figure  4  shows  the  frequency  response  at  angles  of  attack 
for  a  blunt  leading  edge  oscillator  for  a  particular  chamber  and 
inlet  ratio.  The  figure  is  typical  of  the  results  for  each  nozzle 
configuration  and  oscillator  location  and  the  remaining  figures 
similar  to  Figure  4  are  shown  in  Appendix  B.  It  is  important  t<>  note 
on  the  curves  that  the  frequency  response  is  relatively  constant  at 
angles  of  attack  from  -6  to  20  degrees  for  each  run. 

28.  Figure  5  is  a  summary  of  Figures  4  and  Figures  B-l  through 
B-5  of  Appendix  B,  showing  the  oscillator  downstream  location  "d" 
from  the  nozzle  as  parameter  with  frequency  and  inlet  velocity  as 
variables.  The  curves  are  for  Ai/At  =  5.6,  Ac/At  =9.4  and  the 
frequency  values  are  taken  at  near  zero  degree  angle  of  attack.  The 
following  observations  can  be  made  from  this  summary: 

a.  As  the  inlet  velocity  increases  a  linear  divergence  of 
frequency  for  the  family  of  curves  exists.  Within  limits,  the  closer 
the  oscillator  is  to  the  nozzle  exit  the  greater  the  frequency  is 

at  the  higher  inlet  speeds.  The  further  away  the  oscillator  is 
located  the  lower  the  frequency  becomes  at  higher  inlet  air  speeds. 
The  various  downstream  oscillator  locations  indicate  an  optimum 
locaeion  for  constant  frequency. 

b.  The  family  of  curves  can  be  extended  until  a  point  of 
intersection  Is  reached.  Providing  the  linearity  exists  in  this 
region,  the  point  of  intersection  should  indicate  the  threshold 

or  starting  velocity  at  the  systems  natural  frequency.  The  inter¬ 
section  point  is  at  the  calculated  frequency  and  also  indicates 
a  threshold  of  103  knots  inlet  velocity.  More  detail  is  observed 
from  Figure  B-l  of  Appendix  B  showing  a  70  knot  no-go,  95  knot 
imminent  and  lit  knot  go  situation. 

29.  Figure  6  shows  the  same  summary  as  Figure  5  but  the  para¬ 
meter  is  the  inlet  velocity  plotted  against  the  frequency  and 
downstream  oscilla  or  location.  The  point  in  common  for  the  four 
curves  indicates  an  optimum  "d"  location  for  a  constant  frequency. 

For  the  particular  .’Ozzle  configuration  and  Ai/At  Ac/At  ratios,  the 
oscillator  location  \  i  .232  in.  downstream  of  the  nozzle  exit.  The 
intersection  point  of  the  curves  coincides  with  the  calculated 
natural  frequency  (error  of  .4  Hz) . 

30.  Figure  7  shows  the  frequency  response  against  inlet  air 
speed  at  an  oscillator  locution  of  .222  in.  downstream  of  the  nozzle. 
The  frequency  starts  at  the  calculated  natural  frequency  (34.4  H») 
and  is  driven  to  36.0  Hz  at  500  knots  inlet  velocity.  The  oscillator 
location  was  slightly  closer  to  the  nozzle  than  the  optimum  dimension 
(Figure  6)  so  a  slight  driving  of  the  frequency  at  high  inlet  speeds 
is  expected.  The  Ai/At  and  Kc/At  ratios  for  Figure  7  differs 
slightly  from  Figure  5. 
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31.  Figure  7  also  shows  the  starting  velocity  as  105  knots  inlet 
air  speed.  Air  speeds  below  105  knots  did  not  start  the  oscillator 
even  when  rotated  through  angles  of  attack  from  -6  to  20  degrees. 

Mach  1  at  the  throat  for  this  A^/At  ratio  was  calculated  at  117  knots. 

32.  Figure  8  shows  the  frequency  response  with  oscillator  down¬ 
stream  location  as  parameter  for  constant  Ai/At  =  5.1  and  Ac/At  * 

7.40.  The  frequency  starts  at  the  natural  frequency  and  then  diverges 
similar  to  Figure  5 . 

33.  Comparing  the  results  of  Figure  5  to  Figure  8  shows  little 
difference  for_the  frequency  response  for  assumed  constant  Ai/At 
but  different  Ac/At  ratios  at  reasonable  oscillator  locations. 
Placement  of  the  oscillator  too  far  from  the  nozzle  will  eventually 
lead  to  a  non-vibrating  condition  at  high  subsonic  air  speeds.  If 
the  oscillator  is  placed  too  close  to  the  nozzle  a  condition  indicated 
by  Figure  9  is  generated.  The  frequency  response  is  erratic  varying 
from  0  to  55  Hz. 

34.  Figure  10  (summary  of  Figures  B-3,_B-6  and  B-7  of  Appendix 
B)  shows  the  frequency  trends  for  different  Ac/At  ratios.  The 
oscillator  location  was  .280  in.,  near  the  optimum  "d"  dimension  of 
.232.  For  comparison  reasons  only  the  Ai/At  ratioswill  be  considered 
constant  even  though  they  varied  from  5.5_to  4.53.  The  figure  indi¬ 
cates  a  reasonable  frequency  response  at  Ac/At  ratios  of  9.5  to  7.4; 
however,  at  an  Ac/At  of  5.2  the  frequency  has  a  significant  drop-off. 

35.  Tests  were  conducted  comparing  the  response  of  a  rounded 
leading  edge  oscillator  to  a  blunt  leading  edge.  Figure  11  shows 
the  results  for  an  oscillator  location  of  .222"  and  ratios  of  Ai/At  = 
5.07  and  Ac/At  =7.4.  Of  the  two  configurations  the  rounded  oscil¬ 
lator  had  a  30  knot  lower  starting  velocity  than  did  the  blunt  leading 
edge  oscillator.  The  frequency  of  vibrations  at  threshold  velocity 
was  the  systems  natural  frequency  (34.1  H-)  for  both  leading  edge 
configurations.  Both  had  divergence  as  the  inlet  air  speed  increased 
but  more  so  for  the  rounded  configuration.  The  blunt  leading  edge 
went  to  35.7  Hz  and  the  round  leading  edge  to  39  Hz  at  500  knots 
inlet.  (Figure  11  is  a  comparison  of  Figures  7  and  B-8) . 

36.  A  configuration  was  tested  that  had  the  aft  end  of  a  blunt 
leading  edge  oscillator  removed.  A  counterweight  located  out  of 
the  air  stream  was  added  to  dynamically  balance  the  remaining  half. 

The  intent  of  the  setup  was  to  demonstrate  tdiether  or  not  the  FAM 
concept  would  work  with  just  the  front  half  of  the  oscillator  thus 
conserving  space  for  fuzing  applications.  The  blunt  leading  edge 
oscillator  was  i'063  thick  as  compared  to  the  "088  thick  oscillator 

in  the  other,  tests.  A  .203  in.  oscillator  "d"  location  was  used  with 
Ai/At  =4.58  and  Ac/At  =  7.75.  The  results  of  the  test  are  shown  on 
Figure  12.  The  threshold  for  this  configuration  is  estimated  at  250 
knots  with  oscillator  vibration  continuing  to  690  knots  (Mach  1.15) 
inlet  velocity.  Unfortunately  no  real  assessment  of  the  higher 


10 


NOLTR  67-144 


starting  velocity  for  the  one-half  balanced  oscillator  can  be  made. 

Had  a  "063  thick  full  oscillator  been  tested,  a  direct  comparison 
would  be  possible.  The  effect  of  the  oscillator  cross-sectional  area 
on  frequency  response  is  not  known.  The  theoretical  starting  velocity 
was  calculated  at  135  knots  chamber  or  195  knots  inlet.  (Appendix  A.) 
The  starting  vibrations  were  38  H  ,  slightly  lower  than  the  systems 
calculated  natural  frequency  of  40.4  H,.  As  the  inlet  velocity  in¬ 
creased,  the  vibrations  increased  to  42  Hz  for  an  inlet  velocity 
equal  to  Mach  1.  The  chamber  velocity  at  the  leading  edge  was 
calculated  at  400  knots  for  a  free  stream  Mach  1.  After  this  point, 
the  frequency  dropped  to  39  Hz  at  690  knots  inlet  or  475  knots 
chamber  velocity. 

37.  The  rise  time  for  the  ‘J088  thick  blunt  oscillator  to  attain 
its  natural  frequency  was  measured  at  .005  sec.  at  threshold  velocity. 
The  oscillator  "d"  location  was  .280  in.  with  A./A  =  5.6  and  Ac/At  = 
9.5.  The  rise  time  is  for  a  no-load  work  output  condition  and 
measured  from  a  relative  zero  degree  starting  point  to  a  10  degree 
maximum  oscillator  movement  (1/4  cycle) .  The  calculated  rise  time 
for  the  oscillator  to  go  1/4  cycle  at  a  natural  frequency  of 

34.4  Hz  is  .0073  sec.  The  difference  between  the  actual  and 
calculated  rise  times  can  be  attributed  to  1)  the  oscillator  did 
not  start  at  a  0  degree  angle  of  attack,  so  something  less  than  a 
1/4  cycle  was  actually  witnessed  and  2)  the  rise  time  was  measured 
at  tunnel  start-up  where  pressure  surges  of  2  or  3  times  steady  state 
conditions  exist.  It  should  be  noted  that  the  slowest  rise  time 
expected  (no-load)  should  be  based  on  the  natural  frequency 
1/4  cycle  calculation. 

38.  Figure  13  shows  the  chamber  velocity  as  a  function  of  the 
inlet  velocity  with  A^/A^.  as  the  parameter.  The  calculated  threshold 
velocity  is  121  knots  chamber  or  200  knots  inlet  for  Aj/A^-  =  5.6. 

The  actual  threshold  determined  through  tests  was  103  knots  inlet 
(Figure  5)  for  a  blunt  oscillator  at  d  =  .280  in.  The  actual  start¬ 
ing  velocity  almost  coincides  with  Mach  1  at  the  throat.  For  the 
blunt  oscillator  with  d  =  222  and  a  A^/At  =  5.1  (Figure  8)  the  same 

calculated  chamber  velocity  indicates  the  inlet  threshold  should  be 
160  knots.  The  actual  inlet  velocity  was  found  to  be  105  knots 
(Figure  7).  Once  again  the  actual  threshold  approaches  -  the  velocity 
where  Mach  1  is  established  at  the  throat. 

39.  Figure  14  shows  the  typical  change  of  the  static  chamber  air 
density  as  a  function  of  the  chamber  velocity.  The  density  change 
over  the  velocity  range  is  significant  ranging  from  .0022  lb-sec2/ft4 
at  100  knots  to  .0014  lb-sec2/ft4  at  344  knots  chamber  velocity. 

This  density  has  a  bearing  on  the  amount  of  aerodynamic  moment  avail¬ 
able  to  cause  oscillator  vibration  (reference  Equation  8) . 

40.  Figure  15  shows  typical  chamber  velocity  values  at  angles  of 
attack  for  A^/At  =  5.5  and  Ac/At  =  9.5  at  an  oscillator  location  of 

d  =.280  inches.  At  -6  degrees  the  chamber  velocity  started  at 

291  knots,  at  2  degrees  the  minimum  velocity  was  282.4  knots  and  at 

20  degrees  the  chamber  velocity  reached  a  maximum  value  of  302  knots. 
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From  -6  to  +2  degrees  angle  of  attack  the  velocity  relationship  is 
not  linear.  The  velocity  is  influenced  primarily  by  the  blockage 
in  the  tunnel  due  to  the  model  position,  sector  arm,  and  to  a  lesser 
extent  by  a  slight  air  flow  inclination  upstream  at  the  supply  end. 
These  factors  give  a  non-linear,  higher  velocity  reading  at  negative 
and  near  zero  angle  of  attack.  As  the  model  is  rotated  from  +2  to 
positive  angles  of  attack  (up  to>20°)  the  velocity  is  linear.  A 
velocity  value  at  -6°  corresponds  to  a  velocity  value  at  +10°  angle 
of  attack.  At  10°  more  of  the  sector  arm  is  in  the  air  flow  causing 
increased  blockage  that  gives  a  lower  air  speed.  For  summary  figures 
like  Figure  5  the  values  of  chamber  velocity  were  taken  at  near  zero 
degrees  angle  of  attack. 

COMMENTS 

41.  One  theory  explaining  the  flutter  mechanism  that  can  be 
generated  from  the  enclosed  data  will  be  called  "region  of  energy 
gradient. " 


42.  The  basic  premise  for  the  flutter  arming  phenomenon  is  an 
aerodynamic  lift  working  through  a  center  of  pressure  relative  to  a 
pivot  point  and  the  relative  magnitude  of  energy  the  air  possesses. 
The  nozzle  establishes  Sn  energy  region,  i.e. ,  region  of  changing 
velocity  and  pressure,  that  works  on  the  oscillator  through 
aerodynamic  lift.  The  location  of  the  energy  region  is  influenced 
by  the  air  forming  a  nozzle  contour,  aft  of  the  nozzi«v,the  presence 
and  movement  of  the  oscillator  in  the  air  flow,  the  magnitude  of  the 
air  velocity  at  the  inlet,  and  the  magnitude  of  the  back  pressure. 

43.  If  the  oscillator  is  placed  in  a  region  of  undisturbed 
constant  flow  the  aerodynamic  moment  causing  vibration  is  relatively 
low.  This  is  characterized  by  the  placement  of  the  oscillator  too 
far  aft  from  the  nozzle,  out  of  the  energy  region.  The  energy 
gradient  appears  to  move  forward  with  increase  in  velocity  allowing 
even  less  of  the  oscillator  to  be  in  the  energy  region. 

44.  Placing  the  leading  edge  of  the  oscillator  close  to  the 
nozzle  will  place  more  of  the  oscillator  in  the  region  of  energy 
gradient.  A  larger  aerodynamic  moment  is  witnessed  and  in  turn  an 
increasing  frequency  is  established  with  an  increase  in  velocity. 
Midway  between  the  two  extremes  the  frequency  does  not  change  with 
increased  air  velocity.  The  decrease  in  air  density  as  velocity  is 
increased  will  have  a  tempering  effect  tending  toward  a  constant 
frequency  response. 

45.  The  FAM  operation  theory,  "region  of  energy  gradient",  can 
best  be  described  by  an  illustration.  Figure  16  shows  a  qualitative 
diagram  of  the  author's  concept  of  this  theory.  Schlieren 
photography  is  planned  in  the  near  future  to  verify  the  flow  concept. 
Note  that  the  diagrams  show  the  "energy  region"  with  a  stationary 
oscillator  only.  The  flow  pattern  illustrated  is  typical  for  the 
various  physically  defined  Aj/At,  Xc/AT  ratios. 
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a.  Figure  16a  shows  the  relative  location  of  the  energy  region 
and  the  oscillator  location  at  or  near  the  optimum  location.  When 
an  inlet  velocity  (Vj)  is  large  enough  to  establish  Mach  (M)  -  1  at 
the  throat,  M>1  is  established  for  a  small  distance  downstream  of  the 
throat.  Note  that  part  of  the  air  exiting  from  the  nozzle  forms  a 
nozzle  contour  due  to  eddy  flow  and  expansion  conditions.  The  flow 
is  then  decelerated  and  a  second  shock  wave  or  M  -  1  is  established. 
The  region  bounded  by  the  second  shock  wave  to  where  the  flow  has 
reached  steady  state  ( v£ }  is  called  the  "energy  gradient"  region. 

For  a  greater  inlet  velocity  ( V2 )  the  steady  state  region  moves 
forward  (V2'),  accordingly  the  energy  gradient  region  also  moves 
forward,  but  since  the  oscillator  is  at  or  near  the  optimum  location 
it  remains  within  this  region  giving  a  constant  frequency  response. 

b.  Figure  16b  shows  the  oscillator  too  far  aft  from  the  nozzle 
for  effective  use  of  the  "energy  region".  As  the  “energy  region" 
moves  forward  with  increase  of  velocity  the  aerodynamic  moment 
decreases  until  a  non-vibrating  condition  exits. 

c.  Figure  16c  has  the  oscillator  too  close  to  the  nozzle  exit. 
The  second  shock  wave  has  been  altered  somewhat  due  to  the  presence 
of  the  oscillator.  The  "energy  gradient"  exists  over  a  greater 
downstream  distance  allowing  for  a  greater  aerodynamic  moment  and 
increased  frequency  response.  As  the  inlet  velocity  is  increased  the 
shape  of  the  second  shock  wave  becomes  more  concave  and  extends  even 
further  downstream.  The  aerodynamic  moment  increases  to  a  value 
where  the  frequency  response  becomes  very  erratic. 

46.  The  region  of  energy  gradient  for  small  inlet  and  chamber 
differences  (4.53  and  5.2  respectively)  appears  to  be  quite  short. 

Less  of  an  allowable  oscillator  location  deviation  is  permissible. 

At  A^/At  5.6  and  Sc/At  -  9.5  the  oscillator  location  is  not  as 
critical. 

47.  Another  parameter  that  will  influence  the  oscillator 
vibration  is  the  throat  area  relative  to  the  oscillator  cross-section 
area.  Unfortunately  the  data  enclosed  did  not  parametrically  describe 
this  variable.  Although  the  throat  area  did  vary  somewhat  from  test 
to  test  no  conclusions  can  be  drawn. 

CONCLUSIONS 

43.  The  foregoing  test  results  allow  th<=  following  conclusions: 

a.  The  frequency  of  vibration  is  the  natural  frequency  of 
the  oscillator- restoring  member  combination  (no  output  load)". 

b.  Rise  time  from  0  to  the  natural  frequency  for  this 
model  has  an  order  of  magnitude  of  .005  sec  at  threshhold  velocity. 

c.  There  is  an  optimum  oscillator  location  and  nozzle 
configuration  that  gives  constant  frequency  at  all  speeds  above 
threshhold,  but  below  Mach  1  for  free  stream  velocity  and  angles  of 
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ttack.  For  this  particular  model  the  optimized  conditions  are 
hown  in  Figure  17. 

d.  It  should  be  noted  thac  for  a  constant  frequency 
esponse  the  ratio  of  oscillator  "d"  downstream  location  to  the 
iscillator  thickness  is  2.63.  This  ratio  has  been  found  to  hold 
or  other  models. 

e.  Constant  frequencies  at  all  inlet  air  speeds  below 
lach  1  can  be  achieved  even  with  air  densities  varying  by  36%  from 
hat  at  threshhold. 

f.  The  rounded  leading  edge  oscillator  had  a  lower 
hreshhold  velocity  than  the  blunt. 

g.  The  balanced  half  oscillator  will  function  similar  to 
he  full  oscillator  but  the  starting  velocity  is  higher. 

h.  There  are  specific  lower  limits  for  Ac/At  ratios  for 
onstant  frequency  oscillator  response. 

i.  The  method  of  setting  the  restoring  spring  moment  equal 
o  the  aerodynamic  moment  to  determine  the  minimum  velocity  for 
scillator  start-up  is  not  accurate.  It  appears  that  an  inlet 
elocity  that  establishes  Mach  1  at  the  nozzle  throat  is  the 
ritenon  for  this  configuration. 
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FIG.  3  SPACE  LOCATIONS  IN  CONTROL  VOLUME  MODEL 
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VIBRATION  FREQUENCY  (Hz) 
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FIG.  5  FREQUENCY  RESPONSE  AS  A  FUNCTION  OF  OSCILLATOR  LOCATION 
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FIG.  6  OPTIMUM  OSCILLATOR  LOCATION  FOR  SPECIFIC  Aj/At,Ac/At  RATIOS 


CALCULATED  NATURAL  FREQUENCY 
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FIG.  7  FREQUENCY  RESPONSE  AT  NEAR  OPTIMUM  LOCATION 
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FIG.  10  FREQUENCY  RESPONSE  AT  INLET  AIR  SPEEDS  FOR  RELATIVELY  LITTLE  CHANGE  IN 
A  /A  RATIO  BUT  FOR  LARGE  DIFFERENCES  IN  Sc/Af  RATIOS  AT  d=  0.280  IN. 


ROUND  LEADING  EDGE 


MACH  I 
AT  THROAT 


26 


NOLTR  67-144 


NOLTR  67-144 


(b)  OSCILLATOR  LOCATED  TOO  FAR  FROM  NOZZLE 


AIR  FLOW 


<c)  OSCILLATOR  LOCATED  TOO  CLOSE  TO  NOZZLE 


FIG.  16  ENERGY  GRADIENT  CONCEPT 
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CHAMBER  HEIGHT  =  .725  IN 
A. /A  =  .390  =  5.6 
*  .070 

A  /A  =  .700  -  .034  _  „  r 
c  *  - ^070 - 9*5 


A  /\  =  .700 


=  10.0 


NOTE:  .034  CORRECTION  FOR  STATIC  PRESSURE  PROBE 


FIG.  17  OPTIMUM  DIMENSIONS  FOR  CONTROL  VOLUME  MODEL 
FOR  ONE  PARTICULAR  NOZZLE  AND  OSCILLATOR 
CONFIGURATION 


31 


% 


NOLTR  67-144 


APPENDIX  A 


Part  1 . 


From  Reference  (d) ; 

e  =  MtL  - 

R  3  El 


For  the  restoring  member; 


X 


a  —  M«L 
VR  3  El 

but  Mt  =  M/2 

L  =  Z/2 


lM/.2j_AZ/2) 
3  E  I 


but  K  = 


M  _  12  El 

0  Z 
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APPENDIX  A 

Part  2. 

Threshold  velocity  for  one-half,  balanced,  oscillator. 

=  _ K 

cr  rrp  sa  Equation  (8) 

c 

where  S  =  .595  in2  (%  original  oscillator 

area) 


K  =  102.  x  10-2  lb-in 
a  =  .697  in  (3/4  x  .930) 
Pc  =  2.16  x  10~3  lb  -sec2 


228  ft/gec  or  135  knots  chanfoer 


A-2 


VIBRATION  FREQUENCY  (Hz) 
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-1  VIBRATION  FREQUENCY  VS  ANGLE  OF  ATTACK  FOR  THE  CONTROL  VOLUME  WITH 

0.088  IN.  BLUNT  OSCILLATOR  0.188  IN.  DOWNSTREAM  OF  CHAMBER  THROAT 

(A  ),  FOR  INLET  RATIO,  A  /A  -  5.57  AND  CHAMBER  RATIO,  A  /A  =  9.50 
t  \'  t  c  t 


VIBRATION  FREQUENCY  <Hz) 
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FIG.  6-2  VIBRATION  FREQUENCY  VS  ANGLE  OF  ATTACK  FOR  THE  CONTROL  VOLUME  WITH 
0.088  IN.  BLUNT  OSCILLATOR  0.188  IN.  DOWNSTREAM  OF  CHAMBER  THROAT 
(A.),  FOR  INLET  RATIO,  A  /A  -5.79  AND  CHAMBER  RATIO,  A  /A  =  9.94 

t  i'  t  '  C  t 


VIBRATION  FREQUENCY  (Hz) 


FIG .  B-3  VIBRATION  FREQUENCY  VS  ANGLE  OF  ATTACK  FOR  THE  CONTROL  VOLUME  WITH 
0,088  IN.  BLUNT  OSCILLATOR  0.280  IN.  DOWNSTREAM  OF  CHAMBER  THROAT 

(A  ),  FOR  INLET  RATIO,  A. /A  =  5.50  AND  CHAMBER  RATIO,  A  /A  =  9.50 
’  t  c  t 


VIBRATION  FREQUENCY  (Hz) 
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4  VIBRATION  FREQUENCY  VS  ANGLE  OF  ATTACK  FOR  THE  CONTROL  VOLUME  WITH 
0.088  IN.  BLUNT  OSCILLATOR  0 .344  IN  .  DOWNSTREAM  OF  CHAMBER  THROAT 
(A  ),  FOR  INLET  RATIO,  A. /A  =  5.60  AND  CHAMBER  RATIO,  A  /A  =  9.50 


B-4 


VIBRATION  FREQUENCY  (Hz) 


VIBRATION  FREQUENCY  (Hz) 
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FIG.  B-6  VIBRATION  FREQUENCY  VS  ANGLE  OF  ATTACK  FOR  THE  CONTROL  VOLUME  WITH 
0.088  IN.  BLUNT  OSCILLATOR  0.280  IN.  DOWNSTREAM  OF  CHAMBER  THROAT 
(A  ),  FOR  INLET  RATIO,  A. /A  =4.53  AND  CHAMBER  RATIO,  A  /A  =  5.2 
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FIG.  B-7  VIBRATION  FREQUENCY  VS  ANGIE  OF  ATTACK  FOR  THE  CONTROL  VOLUME  WITH 
0.088  IN.  BLUNT  OSCILLATOR  0.280  IN.  DOWNSTREAM  OF  CHAMBER  THROAT 

(A  ),  FOR  INLET  RATIO,  A. /A  =5.17  AND  CHAMBER  RATIO,  A  /A  =  7.40 
*  r  t  c  t 
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